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Numerical results for the transport coefficients of quark gluon plasma are obtained by 
lattice simulations on on 16 3 x 8 lattice with the quench approximation where we apply 
the gauge action proposed by Iwasaki. The bulk viscosity is consistent with zero, and the 
shear viscosity is slightly smaller than the typical hadron masses. They are not far from 
the simple extrapolation on the figure of perturbative calculation in high temperature 
limit down to T ~ T c . The gluon propagator in the confined and deconfined phases are 
also discussed. 

1. Introduction 

In the phenomenological study of quark gluon plasma(QGP), when its bulk properties 
are concerned, the system of quarks and gluons is usually treated as gas or liquid. Then 
the fundamental parameters of QGP such as transport coefficients, are very important in- 
formation. In the high temperature limit, they are calculated based on the perturbation 
theory[[T]-||. However, as temperature T decreases the perturbative calculation breaks 
down. The purpose of this work is to calculate the transport coefficient from the funda- 
mental theory of QCD by lattice simulations in the vicinity of T c . 

The calculation of transport coefficients on the lattice is formulated in the framework 
of linear response theory of Kubof^j4|-|(|. They are expressed by the space time integral of 
retarded Green's function of energy momentum tensors at finite temperature. The shear 
viscosity r\ is expressed as, 
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Similarly the bulk viscosity and heat conductivity are expressed in terms of the retarded 
Green's function of Tn and T41 components of energy momentum tensor. The direct 
calculation of the retarded Green's function at finite temperature is very difficult. Then 
the shortcut is to calculate Matsubara Green's function(C7/3) and then by the analytic 
continuation, we obtain the retarded Green's function at finite temperature. The analytic 
continuation is carried out by the use of the fact that the spectral function of Fourier 
transform of the both Green's functions is the same. For the spectral function we use the 
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following simplest ansatzQ, 

A 7 7 

P(P = 0, W) = — ( ? rj- g - 7 ■ r^- 2) (2) 

71 (m — uj) z + 7^ (m + cjJ 2 + 7^ 

where 7 is related to the imaginary part of self energy and partially represents the effects 
of the interactions. Under this ansatz, the transport coefficients are calculated as, 

o t 27m , , 

where a represents rj, |?7 + C an d X ' T, and a is lattice spacing. 

Notice that if m = or 7 = 0, transport coefficient becomes zero. In order to determine 
these three parameters, at least three independent data points in Gp are required in the 
temperature direction, which means > 6. 



2. Numerical Results on Transport Coefficients of QGP 

It is found that the the fluctuation of Gp is large, and it is a very important problem 
to reduce the fluctuation of Gp. We find that by using the Iwasaki's improved action, 
the fluctuation is much reduced|8|]. Then we apply Iwasaki's Improved action for the 
simulation of SU(3) gauge theory. 

The finite temperature transition point at Nt = 8 with Iwasaki's improved action is 
(3 ~ 2.72 on 16 3 x 8 lattice. We choose our simulation point at f3 = 3.05, 3.2 and 3.3. 
From 0.5 x 10 6 ~ 10 6 Monte Carlo data, we obtain Gp(t) for T n and T 12 . But they have 
still rather large errors. The fit of Gp with parameters by Eq.(2) is done with SALS. The 
shear and bulk viscosities are obtained by these parameters by Eq.(3), and the errors are 
estimated by the Jackknife method. 

The results for shear and bulk viscosities are shown in Fig.l. It is found that the bulk 
viscosity is smaller than shear viscosity and is consistent with zero within errors. This 
is consistent with the result of perturbative calculation at high temperature limit, which 
has been ( = O|0-§. 

We compare our results for shear viscosity with those perturbative formula at high 
temperature limit [|1]-[3|] , r\ = C -T 3 / (—a 2 s loga s ). Where a s is a running coupling constants 
given by, ot s = |y/og(^), A is a scale parameter of QCD. C is a constant depending on 
the method of calculation (0.06 < C < 0.25). By using the relation T = 1/(N T x a), 
perturbative formula is expressed by the lattice spacing a as follows, 

V a 3 = C/(-N 3 a 2 Joga s ) (4) 

The scale parameter A on the lattice is determined by assuming asymptotic scaling 
relation and (3 C at Nt = 8. From the finite size scaling relation reported by the Tsukuba 
group 0, f3 c ~ 2.74 at V = 00 for Nt = 8. And using T c ~ 27QMeV @ for improved 
action, A is determined by, A/T c ~ 1.5. With these quantities the perturbative formula 
Eq.(4) is compared with our calculation as shown in Fig. 2. It is found that our result are 
located not so different from the simple extrapolation of perturbative calculation on the 
figure. But notice that the Eq.(4) breaks down around T/T c ~ 1. 

The shear viscosity in the physical unit are obtained if we assume asymptotic scaling 



Figure 1. Shear and Bulk viscosity 



Figure 2. Shear viscosity from numerical 
and perturbative calculation 



relation for (3 > 2.74 region. The results are shown in Fig.3. They are slightly smaller 
than the ordinally hadron masses. What is the physical effects on the phenomenology of 
quark gluon plasma, when it has shear viscosity with this magnitudes, is a very interesting 
problem. 

In our calculation, Gp of T i4 from which the heat conductivity is calculated, has large 
background and we could not get signal from it. 



3. Gluon Propagators at Finite Temperature 

At first stage of our calculation we planed to calculate transport coefficients both on 
confined and deconfined phase. From the simulation of U(l) gauge theory, it is found that 
the fluctuation of Gp in the confined phase is much larger than those of deconfined phase, 
and we could not determine the Green's functions even with about 1.5 x 10 6 datafTti 



in confined phase. Similar results are obtained in the case of SU(2) gauge theory. To 
find the the difference of Gp between confined and deconfined phase, we study the gluon 
propagators at finite temperature on 8 3 x 4 lattice, because Gp are expressed by the gauge 
invariant combination of gluon propagators. However the gluon propagator itself is gauge 
dependent, and we should fix the gauge. It is done by maximize I = Re(Y^Tr(U fi )) by 



successive gauge transformation fllTf , which is a lattice version of Lorentz gauge. 

The stopping condition for maximizing / is R = \I n+1 — I n \/I n < 10~ 14 . It is found 
that the gluon propagators are very sensitive to R when R > 10~ 9 . With this stop- 
ping condition, we could not find the gauge copies in the deconfined phase. While in the 
confined phase, gauge copies are found. Then in the confined phase, we repeat random 
gauge transformation and gauge fixing processes 5 times, and calculate the gluon propa- 
gators on the gauge fixed configuration with largest values of /. In Fig. 4, the transverse 
gluon propagators are shown for confined and deconfined phase which has a momentum 
of p — (27r/8, 0, 0) and propagate in the ^-direction which is 8 in lattice size. 

In the deconfined phase, the gluon propagator is well represented by the the free Lorentz 
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Figure 3. Shear viscosity in physical unit Figure 4. Gluon propagator in Confined 

and deconfined phase 



gauge gluon propagator on the lattice. While in the confined phase, the propagator shows 



quite different behavior as first found by Nakamura||12|| at the large distance with standard 
action. In Fig. 4 it is found that the effective mass of gluon seems to increase with the 
distance. What is a physical meaning of this behavior should further be investigated, but 
we find that this is a reason why the Matsubara Green's function of energy momentum 
tensor is so noisy in the confined phase. 
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Fig.1 Shear and Bulk Viscosity 
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Fig. 2 Shear Viscosity from Numerical and Perturbative 
Calculation 
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Fig.3 Shear Viscosity in Physical Unit 
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Fig. 4 Gluon Propagator in Confined And Deconfined Phase 



